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INTRODUCTION AND SUMMARY 

This la the tenth quarterly report on Contract NAS 5-27382 
entitled "Spectroradlometrlc Calibration of the Theaatlc Mapper.** 
In this report we describe the results of analyses of TM Images 
acquired on July 8 and October 28, 1984, and of a check of the 
calibration of the 1.22-m Integrating sphere at Santa Barbara 
Research Center (SBRC). 

The results obtained from the In -flight calibration attempts 
disagree with the pre -flight calibrations for bands 2 and 4. Con- 
siderable effort haa been expended In an attempt to explain the 
disagreement. The difficult pc!~v to explain Is that the differ- 
ence between the radiances predicted by the radiative transfer 
code (the code radiances) and the radiances predicted by the pre- 
flight calibration (the pre-flight radiances) fluctuate with spec- 
tral band. Because the spectral quantities we measure at White 
Sands show little change with spectral band, these fluctuations 
are not anticipated. Analyses of other targets at White Sunds 
such as clouds, cloud shadows, and water surfaces tend to support 
the pre-flight and Internal calibrator calibrations. The source of 
the disagreement has not beeu Identified. It could be due to (1) a 
computational error in the data reduction, (2) an incorrect assump- 
tion In the input to the radiative transfer code, or (3) Incorrect 
operation of the field equipment. Items (I) and (2) have been 
checked and recheckad; Item (3) can best be checked by repeat cal- 
ibration attempts, and these are now in progress with measurements 
at Maricopa (ground reflectance about 0.25) on May 20 and at White 
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Sands on May 24. These measurements will be repeated at every 
opportunity during the summer. 

It la of vital laportance that we resolve this problem, not 
slaply to prove or disprove the validity of this in-flight calibra- 
tion technique, but to establish the accuracy and reliability of 
this procedure when It is used for ataospherlc correction. 

Me are pleased to note that Kenneth &. Castle and Carol J. 
Kastner successfully defended their Ph.D. dissertations during this 
reporting period. Me gratefully acknowledge the financial support 
of this contract for their studies. Copies of their dlssertatlors, 
which will be available In late suaaer, will be sent to Dr. John 
L. Barker. 

The relevant data for the October 28, 1984, Land sat V cali- 
bration at White Sands are given in Table 1. 


Table 1. Chuck Site Data, October 28, 1984 


Latitude: 32°55' 

Longitude: 106°22' 

Altitude: 1200 a 
Pressure: 663.7 am 
Temperature: 12.4*C 
Relative humidity: 75 X 
Uae of overpass: 10:09:01 MST 
Solar zenith angle: 52.068° 
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ATMOSPHERE 


The first step In calibrating the dete Is to deteraine the 
ataospharlc components, including the vsrious optical depth values. 
On October 28 the Reagan solar radioaeter (SN 002) was used to 
aake the Langley plot aeasureaents . Measurements were begun at 
6:56 a.a. (all tlaes are MST) at a solar zenith angle of about 84* 
and continued until 12:01 p.a. at a solar zenith angle of 47*. 
Measureaents were aade In nine wavelength regions. These data 
were then run In two Langley plot routines: the Kltt Peak prograa 
and tt'i Reagan prograa. The values froa the two routines agreed 
reasonably well, but the t valut » (see Table 2) used in t.va rest 
of the analysis were those froa the Reagan prograa because it is 
believed to be sore accurate, as It allows different weighting for 
different air aasses. The resulting Langley plots are shown in 
Figure 1. 


Table 2. t ext in Reagan Bands 


text 


Band 

1 (ua) 

Kltt Peak 

Reagan 

2 

0.4025 

0.4888 

0.4996 

3 

0.4417 

0.3654 

0.3710 . 

4 

0.5222 

0.2173 

0.2214 

5 

0.6125 

0.1649 

0.1669 

6 

0.6719 

0.1158 

0.1181 

7 

0.7115 

0.0949 

0.0929 

8 

0.7790 

0.0762 

0.0770 

9 

0.8732 

0.0530 

0.0508 

10 

1.0317 

0.0249 

0.0246 


OMGWAL PAGE' IS 
OP POOR QUALITY 





Figure 1. Reagan radiometer Langley plots for October 28, 1984 
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These total t values were then used to deternlne the Individ- 
ual t components: T ray , and t ocn . The Rayleigh component 
was calculated at the appropriate wavelength ualng the measured 
atmospheric pressure given in Table 1. After the Rayleigh compo- 
nent was subtracted , the Mle component was determined by making a 
quadratic fit to the log t vs. log X data. This yielded the three 
constants a„ a,, and a t in the equation 

log t - a, + a x log X + a, (log x) 1 . 

The values determined vere: a, ■ -1.640, a, - -3.390, a t » -2.935. 
This fit was done with the OZONE program, which also determines 
the amount of ozone by subtracting the calculated Hie and Rayleigh 
components from the total measured optical depths in bands where 
there is ozone absorption. From the ozone absorption coefficient 
at these wavelengths, we calculated the amount of ozone on 
October 28 to have been 182.5 matm-cm. This value was then used 
to calculate T 0 zn the desired wavelengths. The OZONE program 
also gave us the v parameter in the Junge model for the aerosol 
size distribution. The obtained v value of 4.16 seemed a little 
large, so a check was made by running a program that just assumes 
a linear fit to the log t vs. log X data and then v ■ -a, + 2. The 
values for a» and ai thus determined were: a, - -1.55 and a t - 
-2.09. This gave a value of v ■ 4.09. As a check on the sensitiv- 
ity to the value of v, the Herman code was run with v ■ 2.65 and v 
■ 4.09. The resulting values for the normalized radiances at the 
sensor varied less than 5Z. 
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F c*m the parameters determined above, the valuea of the t*!*, 
T ray , and t 0 an componanta in the TM banda vara determined. For 
the TM banda the band llmlte, widths, and cantera uaed vara thoae 
determined by Palmer ualng the momenta method. Theae are aumma- 
rlaed In Tabla 3. In banda 4, 5, and 7 we alao need to account for 
the water vapor and C0 S abaorptlon. Theae were determined by run- 
ning LOWTRAN VI and getting an average water vapor and C0 a 
transmittance over theae banda. The reaultlng tH t O was then 
scaled by tlie measured relative humidity given In Table 1. The 
final t components for each band are summarized in Table 4. These 
were the values used In running the Herman code. 


Table 3. TM Band Parameters 


band 

X| (u«) 

X* (pm) 

AX (pm) 

X c (pm) 

1 

0.4513 

0.5214 

0.0701 

0.4863 

2 

0.5262 

0.6150 

0.0889 

0.5706 

3 

0.6226 

0.6988 

0.0762 

0.6607 

4 

0.7710 

0.9053 

0.1343 

0.8382 

5 

1.564 

1.790 

0.227 

1.677 

7 

2.083 

2.351 

0.268 

2.223 



Table 

4. t Components 



Band 

T mie 

Tray 

T ozn 

TH,0 

TOO, 

1 

0.1360 

0.1420 

0.0047 

0. 

0. 

2 

0.1027 

0.0739 

0.0198 

0. 

0. 

3 

0.0750 

0.0407 

0.0098 

0. 

0. 

4 

0.0401 

0.0156 

0.0011 

0.0454 

0. 

5 

0.0028 

0.0010 

0. 

0.1241 

0.0094 

7 

0.0007 

0.0003 

0. 

0.0805 

0.0035 
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REFLECTANCE 

The next step In the calibration la to determine the sand 
raflactanca In tha six TM bands. These vara determined from 
ernes MMR measurements on the North site (see last quarterly 
report) uelng BaSO* panel 1. This panel had been calibrated by Che 
In April 1984, and it vas recalibrated in band 1 upon our return. 
The change in the calibration vas less than 1 X , so the April cali- 
bration data were used. The reflectance values determined are 
summarised In Table 5. The values used In running the Herman code 
for the six TM reflective bands were 0.425, 0.483, 0.517, 0.559, 
0.351, and 0.129. 

SOLAR IRRADIANCE 

Since the Herman code gives the radiance per unit of exo- 
atmospherlc solar irradiance, we need to determine the solar spec- 
tral lrradiance over the TM bands. These calculations were per- 
formed by two methods, both of which used values for the exo- 
atmospherlc solar spectral Irradiance from Iqbal. In the first 
method, these values were used in the equation 

I E. X U) R(X) dX 


R(X) dX 

■ 0 


where Is the average spectral Irradiance over the band, E,^ is 



8 


Table 5. Absolut* Reflectance 


Use 

Chi 

Ch2 

Ch3 

Ch4 

Ch5 

Cb6 

Ch7 

Be SO. 

9.43 

0.896 

0.891 

0.875 

0.858 

0.816 

0.772 

0.713 

9.53 

0.902 

0.897 

0.881 

0.863 

0.821 

0.774 

0.713 

1st Scan 

9.43 

0.434 

0.487 

0.517 

0.557 

0.512 

0.362 

0.137 

9.44 

0.421 

0.475 

0.507 

0.547 

0.503 

0.351 

0.130 

9.44 

0.427 

0.485 

0.M9 

0.556 

0.510 

0.362 

0.137 

9.45 

0.421 

0.480 

0.514 

0.553 

0.503 

0.348 

0.131 

9.45 

0.427 

0.486 

0.521 

0.564 

0.521 

0.367 

0.134 

9.46 

0.422 

0.481 

0.516 

0.559 

0.511 

0.351 

0.130 

9.46 

0.439 

0.499 

0.533 

0.574 

0.521 

0.362 

0.134 

9.47 

0.438 

0.496 

0.531 

0.575 

0.529 

0.372 

0.141 

9.47 

0.414 

0.474 

0.511 

0.547 

0.498 

0.345 

0.123 

9.48 

0.419 

0.475 

0.507 

0.546 

0.491 

0.332 

0.12e 

9.48 

C.414 

0.468 

0.502 

0.540 

0.481 

0.323 

0.117 

9.49 

0.424 

0.482 

0.517 

0.557 

0.507 

0.351 

0.129 

9.49 

0.415 

0.473 

0.511 

0.559 

0.504 

0.340 

0.125 

9.50 

0.428 

0.490 

0.526 

0.569 

0.525 

0.367 

0.133 

9.50 

0.403 

0.467 

0.506 

0.551 

0.505 

0.343 

0.123 

9.51 

0.434 

0.495 

0.533 

0.581 

0.536 

0.372 

0.136 

Mean 

0.424 

0.482 

0.517 

0.558 

0,510 

0.353 

0.130 

SDEV 

0.010 

0.010 

0.010 

0.012 

0.014 

0.014 

0.006 

BaSO. 

10.01 

0.906 

0.902 

0.886 

0.868 

0.823 

0.775 

0.712 

10.12 

0.912 

0.908 

0.893 

0.873 

0.826 

0.776 

0.710 

2nd Scan 

10.02 

0.427 

0.486 

0.517 

0.553 

0.503 

0.350 

0.130 

10.02 

0.437 

0.496 

0.532 

0.571 

0.518 

0.358 

0.132 

10.03 

0.411 

0.467 

0.501 

0.538 

0.496 

0.348 

0.131 

10.03 

0.420 

0.480 

0.515 

0.554 

0.509 

0.355 

0.133 

10.04 

0.421 

0.476 

0.510 

0.558 

0.514 

0.360 

0.132 

10.04 

0.418 

0.477 

0.514 

0.560 

0.508 

0.347 

0.127 

10.05 

0.436 

0.494 

0.528 

0.571 

0.515 

0.355 

0.129 

10.05 

0.441 

0.498 

0.533 

0.578 

0.532 

0.375 

0.142 

10.06 

0.429 

0.489 

0.527 

0.567 

0.512 

0.354 

0.125 

10.06 

0.422 

0.478 

0.511 

0.549 

0.491 

0.331 

0.122 

10.07 

0.426 

0.482 

0.515 

0.555 

0.494 

0.331 

0.117 

10.07 

0.419 

0.478 

0.513 

0.548 

0.497 

0,345 

0.126 

10.08 

0.407 

0.465 

0.502 

0.550 

0.496 

0.333 

0.121 

10.09 

0.433 

0.493 

0.531 

0.576 

0.535 

0.378 

0.140 

10.09 

0.425 

0.480 

0.514 

0.557 

0.499 

0.337 

0.120 

10.10 

U.430 

0.490 

0.525 

0.569 

0.524 

0.365 

0.131 

Mean 

0.425 

0.483 

0.518 

0.560 

0.509 

0.351 

0.129 

SDEV 

0.009 

0.010 

0.010 

0.011 

0.014 

0.014 

0.007 
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th« Iqtal spectral Irrsdlsnce, and R(X) is tha reaponsivlty from 
"Spectral Characteristics of tha Laudsat Thematic Sensors" by B. L. 
Markham and J. L. Barker (NASA TN-83955). In the second method, 
the solar spectral lrradlance values were mere'.y Integrated over 
the TM band widths as determined by Palmer's moments method. 
Table 6 lists the spectral lrradlance values determined for the 
six bands. Agreement between the two methods was very good, 
differing by less than 0.1Z for five bands and by less than 0.5 Z 
for band 5. 


Table 6. Solar Spectral lrradlance 


Band 

hi < w ■ - 

Method l 

ar” 1 itm"‘) 
Method 2 

Percent 

difference 

1 

1955.475 

1954.182 

-0.066 

2 

1826.889 

1827.331 

0.024 

3 

1544.979 

1543.539 

-0.093 

4 

1042.836 

1043.585 

0.072 

5 

220.186 

219.170 

-0.461 

7 

74.777 

74.788 

0.015 


CODE PREDICTIONS 

The Herman code was run using the t values listed in Table 4 
and the reflectances listed In Table 5. The normal assumptions 
for the particle size distribution were made. The Index of refrac- 
tion was assumed to be n ■ 1.54 - 0.011. The Junge v value was 
4.09. The code was run at the band centers for zenith angles 45* 
and 55*. On October 28 the time of overflight was 10:09:01 a.m., 
giving a solar zenith angle of 52.068*. 


The output radiance was 
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interpolated to this eolar senlth angle and sealed by E^/dS where 
E, x cane iron the "Method 1” column In Table 6, and d ■ 0.9932 AU 
vaa ths eolar distance on October 28. The resulting values are 
the predicted spectral radiances at the entrance pupil and are 
listed in Teble 7. 

Table 7. Herman Code Predictions 


Band 

L x /E, x 

(ar-‘) 

L X 

(W m“* sr“* um“‘) 

1 

0.0784 

135.3130 

2 

0.0842 

155.9754 

3 

0.0931 

145.6780 

4 

0.0927 

98.0695 


MEASURED SPECTRAL RADIANCE 

To compare with the predicted values, vs need the pre-flight 
calibration values for gain and offset and the digital counts for 
our site. ; The digital counts were taken from computer-compatible 
Tape A on October 28, where the site was determined by measuring 
from the road and the helicopter pad. The aver&ge digital count 
over the 16-pixal site was determined for each band. An average 
detector gain and offset (averaged over 16 detectors for the Land- 
sat-5 TM) was determined for each band. This method la appropri- 
ate because there la relatively little variation in the gains and 
offsets over the 16 detectors. These average values were then 
used In the equation 

L X - 


(dc-offset )/galn 
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to determine the aaaaured spectral radiance at the sensor. The 
digital counts In the area of our alts In band 3 are given In 
Figure 2. The average gains and offsets and digital counts used 
and the resulting spectral radiance determined are given in Table 
8 along with the percentage difference froa the predicted spectral 
radiance values in Table 7. 
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Figure 2. Digital counts in band 3. 
Box encloses 16-olxel site. 


Table 8. October 28, 1984, In-Flight Calibration Compared 
with Pre -'light Calibration 


Digital 

counts 

Gain 

Offset 

(W ■'* er“* u* -1 ) 

Z difference, 
code -measured 

measured 

223.250 

15.552 

1.833 

142.37 

f/.l 

171.125 

7.859 

1.690 

215.58 

-27.6 

164.813 

10.203 

1.885 

159.68 

-8.8 

166.375 

10.821 

2.237 

151.69 

-35.3 
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HELICOPTER DATA 

As a verification of the coda, the Caatle redlonetar waa 
flown In tha hallcoptar to obtain aeasuraaents at lntaraedlate 
altitudes. Tha north alta waa aaasurad at 500, 1000, 2000, and 
6000 feat above tha ground level of 3936 feet above sea level. 
The tlae of each aeaaureaant was recorded along with the output 
voltage for each of tha 10 bands. Using the calibration data pro- 
vided by Caatle froa a detector-baaed calibration aada In June 
1984 and updated by Blggar using a source-based calibration In 
April 1985, the digital counts were converted to spectral radiance 
values . 

To coapare these with tne Herman code oredlctlons, we need 
the t coaponents and reflectances for each of the Caatle bands. 
These coaponeots were detaralned in the saae way as for tha TM 
bands, using the paraaeters froa tha Langley plots. Froe the 
LOWTRAN VI run It waa detaralned that CO, la not laportann in any 
of these bands. The aaount of water vapor absorption was deter- 
ained; Castle band 9 is the only one with significant H,0. The t 
coapouents used for the Caatle bands are listed in Tabl' 9. 

We also need the reflectance of the sand lu each of the 
Castle bands. Since this was not aeasured at the site, we used 
sobs sand that had been brought bach previously. It was aaasurad 
in an Integrating sphere reflectoaeter In the Optical Sciences 
Center aeaaureaents lab. Measureaents were aada with varying 
aaounts of water added to atteapt to *laulata the sand coualtlons 
of October 28. These aeasureaents were aade over the wavelength 
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region 0.4 Co 1.0 tim. These were Chen flC Co e sixth-order polyno- 
mial in X, end reflectances were predicted for the four Bernes 
bends within this range and also for the 10 Castle bands for one 
gypsum moisture condition. Since the predicted values for the 
Bernes reflectances agreed with the measured Barnes reflectances, 
the predicted values for the Castle bands were assumed to be 
accurate. The predicted and measured values for the Barnes bands 
and the predicted values for the Castle bands are listed in Table 
10. The Castle band reflectance values were used la the Herman 
code for the intermediate altitudes. Output from the code was in- 
terpolated to the correct solar zenith angle at the time of each 
helicopter measurement, and the exoatmospherlc solar spectral 
lrradlance was determined for each band as above (see Table 9). 
The Interpolated spectral radiance value was then scaled by 
B.^/d*. These values were then compared with the measured values 
from the radiometer. The results ar* given In Table 11. 


Table 9. Castle t Components 


(W m“* sr“ l 


Band 

x c 

T mie 

T ray 

T ozn 

•HjO 

*C0, 

urn” 1 ) 

1 

0.4020 

0.1746 

0.3113 

0 . 

0 . 

0 . 

1555.06 

2 

0.4210 

0.1657 

0.2571 

0.0003 

0 . 

0 . 

1727.96 

3 

0.4403 

0.1567 

0.2136 

0.0005 

0 . 

0 . 

1833.65 

4 

0.5254 

0.1197 

0.1035 

0.0116 

0 . 

0 . 

1893.47 

5 

0.6054 

0.0911 

0.0581 

0.0227 

0 . 

0 . 

1736.33 

6 

0.6621 

0.0746 

0.0404 

0.0096 

0.0016 

0 . 

1540.00 

7 

0.7807 

0.0490 

0.0207 

0.0023 

0.0031 

0 . 

1180.05 

8 

0.8617 

0.0369 

0.0139 

0.0007 

0.0032 

0 . 

987.66 

9 

0.9497 

0.0272 

0.0094 

0 . 

0.4047 

0 . 

782.94 

10 

1.0423 

0.0199 

0.0065 

0 . 

0.0046 

0 . 

679.93 
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Table 10. Beflactancas 


Band 

1 (u») 

Predicted 

Measured 

1 

0.486 

Barnee 

0.453 

0.425 

2 

0.571 

0.508 

0.483 

3 

0.661 

0.537 

0.517 

4 

0.838 

0.550 

0.555 

1 

0.402 

Caetle 

0.444 


2 

0.421 

0.433 

— 

3 

0.440 

0.432 

— 

4 

0.525 


— 

5 

0.605 

0.523 

— 

6 

0.662 

0.537 

— 

7 

0.781 

0.5*7 

— 

8 

0.862 

0.550 

— 

9 

0.950 

0.544 

— 

10 

1.04 

C. 552 

— 


Table 11. Helicopter Data 


Altitude 
(ft AGL) 

Castle 

band 

Zenith 

angle 

Radiance (W ■'* sr“‘ u* _ ‘) 
Measured Predicted 

Percent 

difference 

500 

2 

53.633 

103.398 

115.788 

-10.70 


3 

53.621 

114.510 

125.423 

-8.70 


4 

53.610 

152.852 

151.819 

0.68 


5 

53.597 

160.079 

154.450 

3.64 


6 

53.584 

152.993 

146.267 

4.60 


7 

53.567 

134.337 

117.214 

14.61 


8* 

53.555 

112.457 

99.875 

12.60 


9* 

53.542 

39.123 

30.088 

30.03 


10* 

53.531 

63.444 

69.812 

-9.12 

1000 

2 

53.174 

106.344 

117.709 

-9.66 


3 

53.161 

118.445 

127.337 

-6.98 


4 

53.148 

156.382 

153.746 

-1.71 


5 

53.137 

160.698 

156.303 

2.81 


6 

53.124 

150.612 

147.945 

1.80 


7 

53.107 

132.112 

118.496 

11.49 


8* 

53.095 

113.076 

100.950 

12.01 


9* 

53.083 

38.643 

29.605 

30.53 


10* 

53.071 

64.828 

70.539 

-8.10 
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2000 


6000 


6000 


2 

49.010 

119.243 

132.723 

-10.16 

3 

49.003 

130.054 

143.077 

-9.10 

4 

48.994 

171.030 

170.949 

0.04 

5 

48.986 

176.328 

173. J3 

1.86 

6 

48.979 

166.679 

163. ±00 

2.19 

7 

48.970 

142.726 

130.064 

9.74 

8* 

48.963 

124.532 

110.625 

12.57 

9* 

48.954 

41.855 

32.819 

27.53 

10* 

48.946 

67.320 

77.155 

-12.75 

2 

51.315 

114.425 

128.226 

-10.76 

3 

51.306 

124.252 

137.722 

-9.78 

4 

51.295 

160.373 

162.493 

-1.30 

5 

51.284 

171.306 

163.864 

4.54 

6 

51.274 

156.897 

154.544 

1.52 

7 

51.264 

133.824 

123.194 

8.63 

8* 

51.252 

113.802 

104.866 

8.52 

9* 

51.243 

32.626 

24.752 

31.81 

10* 

51.232 

66.723 

73.197 

-8.84 

2 

51.188 

112.452 

128.660 

-12.60 

3 

51.177 

124.252 

138.183 

-10.08 

4 

51.166 

163.343 

163.017 

0.20 

5 

51.156 

171.978 

164.380 

4.62 

6 

51.145 

155.029 

155.013 

0.01 

7 

51.134 

138.675 

123.560 

12.23 

8* 

51.125 

114.595 

105.167 

8.96 

9* 

51.114 

32.626 

24.882 

31.12 

10* 

51.103 

65.089 

73.406 

-11.33 


NOTES: Values for the last three bands (*) are questionable 
because the Instrument beater was not working. The first two 
bands are questionable because of their proximity to strong 
Fraunhofer lines. Two runs were made at 600C ft. 


DIFFUSE-TO-DIRECT MEASUREMENTS 

As a second verification of the code, the ratio of diffuse to 
direct Irradlance was measured with the Barnes radiometer. This 
was done by using the Barnes radiometer to view a BaS0„ panel that 
was alternately shaded and Illuminated directly by the sun. The 
shadowing was done by a large styrofoam panel covered with black 
felt on top of a 3.7-m-long pole. The styrofoam panel blocked 
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•bout 0.028 «r or 0.4 X of the sky. The measurement of the shaded 
panel gave the diffuse path radiance, and the difference between 
the shaded and unshaded panel gave the direct solar ir radiance. 
The ratio of these two aeasureaents was then coapared with the 
predictions froa the radiative transfer cude. The aeasureaents 
were aade throughout the morning except during the tiae the 
Barnes was being used for sand reflectance aeasureaents. Thus, 
they covered a wide range of solar zenith angles. Figure 3 sum- 
aarizes the data in graphical fora for the four TM reflective 
bands as a function of zenith angle. Except for the values around 
zenith angle 55°, the pluses represent measured data and the 
squares represent Herman code predictlcns, with the longer wave- 
length bands having the lower values for a given zenith angle. 
These values were all calculated assuming a refractive index of 
1.54 - O.Gi_. The measured band 1 values fall below the predicted 
values by approximately 15 X. The agreement in bands 2 and 3 is 
fairly good, with measured values slightly exceeding the predicted 
values. In band 4 the measured values are too high by about 35 X. 
Note that these differences are in the same direction as the 
differences between the helicopter measured and predicted radiance 
values for similar wavelengths. 

To test the sensitivity of the dlffuse-to-direct retio to the 
assumed refractive index, we varied both the real and imaginary 
parts separately and looked at the ratio. These values are plot- 
ted close to 55° (all were calculated at 55°, but the ordinate was 
displaced slightly so the various points would not overlap). The 
pluses represent variations in the imaginary component with the 


PREDICTED 



Figure 3. Dif f use-to-direct ratio as a function of solar zenith angle. 
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real part bald at 1.54. The laa ginary value* are: 0.0. -0.005, 
-0.015, and -0.02 tor each of the four band*. A* expected, In each 
caae the run with no abeorptlon (1 ■ 0.0) ha* tha highest ratio. 
Tha dot* represent variation* in tha real refractive Index with the 
Imaginary part held at 1 ■ -0.01. The value* for n were: 1.5, 
1.54, 1.75, and 2.0. The ratio 1* higher for the higher n values. 
I* can be seen that the ratio of diffuse to direct Is sore sensi- 
tive to absorption than to changes in the real part of the Index. 
As can be seen on the plot for each band (except posalbly band 4), 
by varying the Index within these Halts we could force the pre- 
dicted ratio to aatch the measured one, but for each band the 
values would be different. 

The Herman code runs with varying Index were also examined to 
see the effect on the radiance at the sensor. Only in band 1 
could the pre-flight radiances be matched to the code output. For 
the other bands, either values of 1 > 0 (If n - 1.54) or a real re- 
fractive Index greater than 2 (for 1 - -0.01) was required. In 

band 1 a aatch could be made with n - 1.54 and 0 > 1 > -0.005 or 

with 1 - -0.005 and 1.75 < n < 2.0. 

AUXILIARY DATA 

In an attempt to understand the above caHbratlon results, we 
compared the data from July 8 and October 28. In addition to the 
data for the sand areas, digital count values were obtained for 
water and cloud and cloud shadow regions near the White Sands 
site. The water and cloud shadow values were used to compare to 
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path radiance predictions fro* the Herman code. Table 12 contains 
these data for the four TM reflective bands (band 1 was saturated 
on July 8) as well as for ratios between the bands. Table 13 con- 
tains band-by-band ratios of the October 28 values to the July 8 
values. The next section explains the rationale behind some of 
the comparisons and some of the results. 

Variations due to the atmouphere are expected to cause only 
up to around 10 Z variations in the radiance at the sensor. Thus 
the values for the two days as predicted by the atmospheric code 
should not vary by more than a few percent. However, the abso- 
lute values are different owing to the different solar zenith 
angles at the time of overpass for these days. For this reason, 
we looked at the ratios between the bands, to cancel out the 
cosine dependence on solar zenith angle. 

As seen in Tables 12 and 13, the band ratios of the code 
predictions agree very well for the two days. The largest varia- 
tions are in band 4 and are directly attributable to the different 
amounts of water vapor absorption on the two days. In addition to 
listing the code predictions, Table 12 lists the digital counts for 
the sand, water, cloud, and cloud shadow areas as veil as these 
digital counts converted to radiance values by means of the pre- 
flight calibration offsets and gains. Of particular interest are 
the water values in band 4. Owing to the large absorptance of 
water in this band, this value is expected to give a comparison to 
the path radiance as computed from the code for a zero reflec- 
tance region. On October 28, the water value was 9.0 W sr“ l . 
The value obtained from the code for a reflectance path radiance 
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Table 12. Code Predictions, Pre-Flight Radiances, Internal Cali- 
brator Radiances, and Digital Counts, July and October 1984 


Band 


1 

2 

3 

4 

2/3 

3/4 

2/4 




July 8, 1984 





• 


Code 

Prediction 




Total 

radiance 

— 

266.3 

248.2 

158.3 

1.073 

1.568 

1.682 

from 

sand site 








Path radiance 

— » 

37.68 

24.8 

10.5 

1.516 

2.362 

3.581 




Pre -Flight Radiance 




Sand 


— 

251.3 

228.4 

180.7 

1.100 

1.26-4 

1.391 

Water 


— 

57.6 

20.7 

10.0 

2.785 

2.080 

5.794 

Cloud 


— 

261.2 

228.5 

184.6 

1.143 

1.238 

1.415 

Cloud 

shadow 

— 

18.8 

10.9 

16.3 

1.725 

0.669 

1.153 



Internal Calibrator 

Radiance 



Sand 


— 

270.1 

241.3 

185.4 

1.119 

1.302 

1.457 




DLgltal Counts 




Sand 


— 

199.2 

234.9 

197.8 

0.848 

1.188 

1.007 

Water 


— 

47 

23 

13 

2.044 

1.770 

3.615 

Cloud 


— 

207 

235 

202 

0.881 

1.163 

1.025 

Cloud 

shadow 

— 

16.5 

13.0 

19.9 

1.269 

0.653 

0.829 




October 28, 1984 







Code 

Prediction 




Total 

radiance 

155.3 

156.0 

145.7 

98.1 

1.070 

1.485 

1.590 

from 

sand site 








Path radiance 

44.4 

28.9 

16.6 

4.9 

1.741 

3.388 

5.898 




Pre -Flight Radiance 




Sand 


142.4 

215.5 

159.7 

151.7 

1.349 

1.053 

1.421 

Water 


50.3 

30.9 

37.4 

9.0 

0.826 

4.156 

3.433 

Cloud 

1 

Sat. 

249.1 

189.3 

168.0 

1.316 

1.127 

1.483 

Cloud 

shadow 1 

27.2 

30.7 

16.8 

13.6 

1.827 

1.235 

2.257 



Internal Calibrator 

Radiance 



Sand 


155.5 

232.5 

170.1 

157.3 

1.367 

1.08i 

1.478 




DLgltal Counts 




Sand 


223.3 

171.1 

164.8 

166.4 

1.038 

0.990 

1.028 

Water 


80 

26 

40 

12 

0.650 

3.333 

2.167 

Cloud 

1 

Sat. 

197.5 

195 

184 

1.013 

1.060 

1.073 

Cloud 

2 

204 

159.6 

163.5 

169.8 

0.976 

0.963 

0.940 

Cloud 

shadow 1 

44.1 

25.8 

19.0 

17 

1.358 

1.118 

1.518 

Cloud 

shadow 2 

46 

25 

19 

17 

1.316 

1.118 

1.471 
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Table 13. Ratios of October 28 to July 8 

Measureaeuts 


Band 

2 

3 

4 


Code Predictions 


Total i 

radiance 0.586 

0.587 

0.620 

froa 

sand alts 




Pr a -Plight Radiances 


Sand 

0.858 

0.699 

0.840 

Water 

0.536 

1.807 

0.905 

Cloud 

0.954 

0.828 

0.910 

Shadow 

1.633 

1.051 

0.834 


Digital Counts 


Sand 

0.859 

0.702 

0.921 

Water 

0.553 

1.739 

0.923 

Cloud 

0.954 

0.830 

0.911 

Shadow 

1.564 

1.462 

0.854 


Internal Calibrator Radiances 

Sand 

0.861 

0.705 

0.848 


of zero was 2.24 W a -1 sr“‘. An atteapt was aada to force the 
code value to agree with the water value by increasing the 
scattering optical depth coaponents. Both the and T ray com- 

ponents needed to be increased by a factor of 4 to get the path 
radiance to agree with the aeasured value. A check was aade to 
see how auch effect this would have on the radiance at the sensor 
for the sand reflectance. As expected, the sensor radiance value 
decreased, but only by a few percent. Note that the decrease is 
opposite in direction to the increase of 35 X as required by the 
calibration result. Thus, even by increasing the path radiance 
values to natch the aeasured ones (deterained froa the ’’zero'* 
reflectance region) we cannot explain the difference between the 
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code result and tha pre-f light calibration result, and In fact the 
change la opposite to the one required. 

The cloud shadow regions were examined because many people 
attempt to use them to determine path radiance values. Comparison 
shows that the cloud shadow values In band 4 are even higher than 
the water values, which were much higher than the code path radi- 
ances. This Indicates that cloud shadows nay not be a good way to 
determine patn radiance values. It Is Interesting to note that for 
October 28 two cloud shadow regions were examined and the numbers 
for the two are very close. 

Also included in Table 12 are Internal calibrator updates to 
the sand radiance values. Barker has Indicated that comparisons 
between pre-flight calibration-based radiance values and their 
updated values using Internal calibrator data may be unreliable. 
However, the internal calibrator Is believed to be a good Indicator 
of in-flight focal plane and electronics changes. In this respect 
it can be used to compare how much the response of the TM bands 
has changed from July 8 to October 28. 

Other ratios are presented In Table 12, for example for cloud 
regions. It Is known that cloud reflectances are fairly flat over 
the visible region. Thus, any differences In the cloud radiance 
should be due to Rayleigh and ozone scattering above the dol'd 
height and to the amount of solar loading. However, to avoid 
saturation, we had to select regions neer the edges of the clouds, 
and this could Invalidate our assumption of nearly spectrally flat 
reflectance. These and the rest of the values In the table ere 
presented without further Interpretation. 
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INSTRUMENT RECALIBRATION 

Another investigation that has bssn performed is tha recall- 
bratlon of tha Integrating sphere at SBRC that vis used to cali- 
brate tha TM sensors. In August 1984, Fulton at SBRC recalibrated 
the Integrating sphere by using the sane setup as originally used. 
This setup consisted of a standard laap Illuminating a Halon panel 
that was imaged onto the entrance silt of a monochromator with a 
photomultiplier tube following the exit slit. A rotatable mirror 
In the system allowed the Integrating sphere to alternately Irra- 
diate the aouochroaator'a entrance slit. Thus the radiance from 
the Integrating sphere relative to that of the Halon panel was 
determined wavelength by wavelength for wavelength Intervals of 
roughly 50 nm. In March 1985, the Castle radiometer (097), origi- 
nally calibrated with reference to a QED detector and an electri- 
cally calibrated pyroelectric radiometer (ECPR) In June 1984 at the 
Optical Sdeuces Center, was recalibrated at SBRC In bands 3 
through 8 using the same lamp and Halon panel. Finally, the In- 
tegrating sphere was calibrated with the Castle radiometer by com- 
paring the readings from the standard lamp and Halon panel with 
those from the Integrating sphere. In addition, the Castle radiom- 
eter was recalibrated by Bigger In April 1985 using a standard 
lamp and a BaSO, panel at the Optical Sciences Center. The re- 

sults of these recallbratlons and comparisons are listed In Table 
14. It can be seen that the two calibrations of the Integrating 
sphere differ by less than 5Z, the smallest differences being for 
the longer wavelengths. 


f 


» 
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Table 14. Castle (097) and Intagrating Sphara Calibration 


Kasponsivity Intagrating Sphara Fadianca 

(V cn 1 ar ■W~ ,> ) (aW ca"* ar“Q 


Band 

Can tie Bigger Vltaan 
(QdD) (BaSOJ (Halon) 

Castle, 

Lg** 

Monochroaator 
L. (L b -L c )/L b (X) 

3 

17.14 

15.214 

16.885 

0.0 J02 

0.0315 

4.1 

4 

25.06 

24.369 

24.575 

0.0819 

0.0844 

2.9 

b 

36.50 

35.421 

35.256 

0.1506 

0.1545 

2.5 

6 

40.76 

39.943 

39.655 

0.1846 

0.1886 

2.1 

7 

55.73 

55.088 

54.345 

0.2760 

0.2776 

0.6 

8 

45.86* 

45.441 

44.345 

0.3685 

0.3650 

-1.0 


*Thia calibration vaa dona with an ECPR. 

**These valuas ara froa the Hltnan :allbratlon (preceding colunn). 


The a a calibration procedures and results will be described In 
■ore detail by Sandra Wltaan in her M.S. thesis, which should be 
completed this fall. The sphere recalibration work was supported 
by grant NAG5-196. We wish to express our appreciation to Janes 
Young and Linda Fulton of SBRC for their enthusiastic cooperstion 
in this work. 

CONCLUSIONS 

Several laportant conclusions can be drawn from the results 
reported here. 

(1) We have found no grounds to question the validity of the 
pre-flight Integrating sphere absolute calibration. In fact, ig- 
noring any degradation with tlae of the output of the sphere, the 
pre-flight calibration appears to agiea with our callbxatlon to 
within 5Z. 

(2) The internal calibrator data in band 2 indicates a change 
of 7.3 X and 6.9 X froa pre-flight to October 28, 1984, and to July 
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8, 1984, respectively. The chsnge fro* July to October Is there- 
fore only 0.4 X. 

(3) The helicopter dsts for Ceetle bends 3 through 7 (0.44 to 
0.78 mb) et en altitude of *''00 ft ACL show differences between 
code-predicted end aeasured radiances of ft on -10 X et 0.44 mb to 
+12X at 0.78 mb. The direction of this trend is in agreeaent with 
the differences between the code and aeasured results for July 8 
and October 28, 1984. The difference between the radiances at 
6000 ft ACL and space is predicted to be less than IX for TM 
bends 3 end 4 for the cond'tlons at White Sands. 

(4) The above results ere at odds with the code-predicted 
results for the TM calibration of October 28. The worst case TM 
band 4 difference of 35 X could be explained by a failure of the 
detector heater in the solar radioaeter. However, this would not 
explain the discrepancy in TM band 2. 

(5) There is raaarkably good agreenent between the TM bend 
code-predicted ratio values (see Table 12) for the July 8 and 
October 28 calibrations. It is difficult to explain how this can 
occur while at the sane tine the TM band 4 value for July 8 
differs froa the pre-flight vslue by only 12 X and the October 28 
calibration Indicates a difference of 35 X. 

(6) According to our results, the use of water in TM band 4 
or cloud shadows as sero radiance areas can lead to large errors 
in estivating path radiance. If this is verified by later observa- 
tions, it will have serious laplicatlons for ataospher'.c correction 
aethods that depend on this technique for deteraining path 


radiance 



